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Abstract. We give an overview over recent studies on the model of Active Brownian Motion (ABM) coupled 
to reservoirs providing free energy which may be converted into kinetic energy of motion. First, we present 
an introduction to a general concept of active Brownian particles which are capable to take up energy 
from the source and transform part of it in order to perform various activities. In the second part of our 
presentation we consider applications of ABM to ratchet systems with different forms of differentiable 
potentials. Both analytical and numerical evaluations are discussed for three cases of sinusoidal, staircase- 
like and Mateos ratchet potentials, also with the additional loads modeled by tilted potential structure. 
In addition, stochastic character of the kinetics is investigated by considering perturbation by Gaussian 
white noise which is shown to be responsible for driving the directionality of the asymptotic flux in the 
ratchet. This stochastically driven directionality effect is visualized as a strong nonmonotonic dependence 
of the statistics of the right versus left trajectories of motion leading to a net current of particles. Possible 
applications of the ratchet systems to molecular motors are also briefly discussed. 
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PACS. 05.40.-a Fluctuation phenomena, random processes, noise, Brownian motion 
dynamics and chaos 



05. 45. -a Nonlinear 



1 Introduction 

The study of mechanical systems with support of energy 
goes back to the investigations of Helmholtz and Rayleigh 
on the origin of sustained oscillations. Generalizations of 
the "active friction" introduced by Rayleigh found many 
applications including the concept of the Active Brown- 
ian Motion (ABM) which extends the notion of standard 
Brownian motion as studied by Einstein, Smoluchwski, 
Fokker, Planck and others [1] to the field of driven mo- 
tions [lllSllllElH] including a developing theory of swarm- 
ing motions [7ll8l[9lfTT1fT2]. 

In this work we first introduce a model of Langevin dy- 
namics coupled to energy depot dynamics and study basic 
properties. The inclusion of the depot dynamics should 
model the general observation that external energy which 
is needed for acceleration of motion is in most cases stored 
first in a depot or reservoir (a kind of an energy tank) and, 
only in a second step, becomes converted by a kind of mo- 
tor into motion. This storage and subsequent conversion 
of the energy into mechanical work is modeled here in 
the simplest possible way by a balance equation. Next we 
study more specific applications to transport problems on 
Hamiltonian ratchets and discuss applications to molecu- 
lar motors. 



Since the fundamental work of Marian Smoluchowski 
[T] the problem of transport on ratchets is under a con- 
stant debate. Some of the most interesting applications 
are related to biological problems (see e.g. Jiihcher and 
Prost. [T3lll4j ). Here we will discuss several problems re- 
lated to the ABM of particles on ratchet potentials. Fur- 
ther we will discuss possible use of ABM in modelling 
the functioning of ATP and ADP in cells and the related 
transport mechanisms, having in mind possible applica- 
tions to biological systems as e.g. proton pumps and elec- 
tron pumps. Some further relations with the stepping mo- 
tor described e.g. by Bier P^^fTf] are also proposed. 



2 Model of Brownian motion coupled to 
energy reservoirs 

2.1 Coupling between Langevin dynamics and energy 
depot dynamics 

We postulate dynamics of Brownian particles as deter- 
mined by the Langevin equation according to a model 
proposed by Schweitzer et al. [4J: 
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r, =v,; mv, + VC/(r) = F, + V2^e(t) (1) 



2 



Alessaiidro Fiasconaro et al.; Active Brownian Motion Models and Applications to Ratchets 



where U{r) is the potential of the conservative forces and 
\/2D£^{t) is a stochastic force with strength D and a S- 
correlated time dependence: 



(2) 



The driving forces on the RHS of Eq.((T]) are expressed in 
the form |4j 

Fi = — m7v + mde{t)v (3) 

where the first term stands for the dissipative force and 
the second is responsible for the acceleration of movement 
due to the conversion of the depot energy into kinetic en- 
ergy of motion. Under the equilibrium conditions, the fric- 
tion coefficient 7 (here defined as the velocity dependent 
function 7(v)) is a constant directly related to the noise 
strength D by the Einstein relation D = uikBT^. Depot 
energy dissipation and coupling of energy reservoir to the 
kinetic degrees of freedom give rise to the time-dependent 
balance equation for e{t) 



e{t) ^ q - ce{t) - de{t)v^ 



(4) 



Here q is the take-up of energy term and ce(t) describes the 
internal dissipation in the reservoir which is assumed to 
be proportional to the depot energy e(t). The conversion 
of depot into kinetic energy of motion is controlled by the 
rate d and depends (quadratically) on the actual velocity 
of the particle. The overall time variation of the mechan- 
ical energy of the particle can be derived from Eqs.(l-4) 
and yields 

= -^me{t) + m{q - ce{t)) - rnvfj + V2D^{t)vi (5) 

In an adiabatic approximation (assumed rapid relax- 
ation of variations in e{t) to its stationary value, i.e. e{t) = 
0) we may substitute the energy e{t) in Eq. ^ by its sta- 
tionary value and get a general form for the dissipative 
force 

Fi = -m-f{vf)wi (6) 

The function 7(u^) denotes a velocity-dependent friction, 
which in our model has a negative part. Accordingly, the 
depot model [4l[6] for the energy supply leads to 



7(v^) 



7- 



dq 



dv"^ 



(7) 



where c, d, q are certain positive constants characterizing 
the energy flows from the depot to the particles. Depend- 
ing on the parameters 7, c, d and q the dissipative force 
function may have one zero at v = or two more zeros at 



when the bifurcation parameter 

7 



(8) 



(9) 



becomes positive. In this case a finite characteristic ve- 
locity uq exists which determines an attractor of motion 
and one speaks of an active Brownian-particles motion. 
For |v| < Wo, the dissipative force is positive, i.e. the par- 
ticle is provided with additional free energy. Hence, slow 
particles are accelerated, while the motion of fast particles 
becomes damped. Note that in the case of thermal equi- 
librium systems we have 7(v^) = 70 = const. The adia- 
batic treatment of the equations proposed by Schweitzer 
et al. ^ found many applications to problems of active 
Brownian motions as e.g. swarm dynamics 19,21,22.i25]. 
However, in some cases, and in particular for applications 
to ratchet problems the adiabatic approximation under- 
stood as mere adiabatic elimination of fast variables from 
the dynamic equations of motion may lead to improper 
conclusions as will be discussed in the forthcoming para- 
graphs. 



2.2 Free Brownian particles and the action of a 
constant external force 

Let us first discuss the free motion of active particles in 
a one or two-dimensional space, r — {xi,X2}. Under the 
condition of a quasi-stationary depot, e{t) — 0, the long 
time (stationary) probability distribution function of the 
Fokker-Planck equation 



dP(r,v,t) ^^dP{r,v,t) ^ VU(y) dP{r,v,t) 



dt 



dr 



d_ 
dv 



7(r, v)vP{r, v, t) + Dy 



m dv 
dP{r,v,t) 



dv 



(10) 



corresponding to the Langevin equation Eq.(l) with U (r) 
can be easily found and reads [B] 



dv^ 



2\ (9/2£'^.) 



exp 



7 2 

V 



2D 



(11) 



where = D/w? (cf. Eqs. (1,2)) stands for the diffusion 
coefficient in the velocity space. Accordingly, the mean 
square displacement 



<(ri(t)-ri(0))2 >= 

dti / dt2 < v{ti)v{t2) > 
Jo 



(12) 



can be evaluated and in the limit of a strong, supercritical 
influx of energy (q 00) and for times much longer than 
VQ/2Dy leads to an approximate expression 



< (ri(t) - ri(0))2 -JL(2i) = 2D,fft : 



4^2 



2vQm 
D 



-t. 



(13) 



So far the usual Brownian motion and the active Brow- 
nian motion seem to behave in a quite similar way: we 
see that the dispersion of the displacement grows linearly 
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with time resembling typical character of the diffusive mo- 
tion with Deff — v^/Dy . However, the prefactor 2Z?e// 
is completely different from the value 2kBT{rn'y)~^ which 
rules the relation < {ri{t) - 7-1(0))^ >= 2kBT{m-f)~'^ for 
a standard Brownian diffusion. In case of ABM, due to the 
inverse-proportional dependence on the noise intensity D 
(cf. Eq. (1)), weak noise gives rise to large mean square 
displacement which is one of the peculiar properties of the 
motion [3]. 



The solutions may be found graphically (see Fig.[T]). If we 
assume that dissipation c and the slope a are so small that 
ac — —Fqc may be neglected, the solution reads 

Altogether, in the general case we have a cubic equation 



From now on we will consider only one-dimensional 
problems. First we analyze the deterministic equations for 
several instructive special cases (for simplicity of deriva- 
tions we assume m = 1): 

x{t) + U'{x) + [7 - de{t)]i{t) = Fo (14) 

e{t) ~ q + ce{t) + de{t)x^ {t) ^ Q (15) 

Here Fq stands for a possible tilt of the potential. As a 



Constant Force Solution 



zeroth-order approximation we first neglect in Eq. (|14[) the 
term U'{x) = and we adiabatically eliminate the energy 
term from the second equation by approximating e = 0: 



e{t) = ^ e 



X — V 

[7 - de{t)]v + Fo 
g 

: -I- dv^ 



(16) 



When treating the set of the above equations by use of the 
equilibrium condition (vanishing of the dissipation term in 
the second equation) which requires eg — j/d we assume 
that the stability flux of energy counterbalances the fric- 
tion term: 




Fig. 1. Stationary solution: The force F{v) — Fo — —a as a 
function of the velocity for 7 = 0.2 and 2 fixed d- values (d = 0.4 
and d = 0.1 with c = 0.1) and the asymptotic curve with c — 0. 
The comparison is made with the Stokes law F{v) — 71;. The 
solution to Eq. (I19p representing free motion corresponds to 
the roots on the axis F{v) = 



X = V 

v = Fo 

_ 7 _ q 
^0 - - .^^,,2 



c + dvQ 

The above set of equation implies now that Fq = 0. 

We will take later this exact solution (with Fq = and 
U'{x) = 0) as the starting point of a perturbation theory. 
The bifurcation parameter of our problem is (3. For /3 > 0, 
the system is driven to non-equilibrium states (note that 
this case requires pumping of energy from the reservoir) 
and has all together three stationary states of the velocity: 
V = V = vq and v — — wq- 

We consider now the case of a constant tilt by means 
of the additive force 

Fo = -a. (18) 

Here a is the slope of an equivalent potential Uo{x) = ax. 
We will mostly focus on a positive slope a > 0, Fq < 0. 
Let us first consider the case in the absence of an addi- 
tional ratchet potential i.e. U{x) = 0. This problem still 
admits an exact solution. Without an energy flux q from 
the reservoir, the particle would fall down (if a > from 
right to left) . Including the reservoir provides the possibil- 
ity of uphill motions. The condition of stationary motion 
under the action of this force leads to the cubic equation 

7duo ^ Fodv'o + {cj — qd)vo — Fqc ~ 



for the stationary velocities, which is easy to solve nu- 
(17) merically, and in some special cases, also analytically. Fol- 
lowing the structure of the linear potential Uq, the down- 
hill motion exists in all cases and for our standard choice 
a > 0; Fo < the downhill motion is directed to the left. 
Remarkably, in the case of positive energy input g > also 
a stationary uphill motion may exist vq > 0, provided the 
force driving downhill is not too large. For example, if 
a = l;Fo = —1, the trivial downhill solution is wo — —5 
and the stable uphill solution is vq = +0.5 for d = 1; for 
d = 0.3 no uphill solution exists. It is interesting to men- 
tion that the uphill motion may exist even without any 
ratchet effects provided the driving is sufficiently strong. 
Still the question remains, what is the influence of the 
ratchet potential on the directionality of transport and 
what is the mechanical efficiency of the system. 
This brief examination allows us to conclude about the 
minimum set-up conditions for construction of a ratchet- 
type potential in which the unidirectional current can be 
obtained: 

(i) the average value of the flatter slope of the potential 
should be in the range where the uphill motion is possible, 
and 

(ii) the average value of the steeper slope should not allow 
the uphill motion. 

(19) Under these conditions the particle can go uphill from left 



4 



Alessandro Fiasconaro et al.; Active Brownian Motion Models and Applications to Ratchets 



to right and, since the motion backwards is not possible, 
we get a unidirectional movement. We will use this con- 
struction as a rule of thumb in order to find the conditions 
for directed transport in one dimensional periodic struc- 
tures. 



2.3 Effect of white noise in the mechanical equations 
with constant force 

Having in mind biological applications where typically 
some transfer from chemical to mechanic or electric en- 
ergy appears in the presence of random fluctuations, we 
will study now ratchets which are connected to an energy 
reservoir under the influence of noise. As a generalization 
of Eq.dllD 



Sinusoidal Potential 




0.5 1 1.5 2 2.5 3 



X 

Fig. 3. Sinusoidal potential with parameter h — 1.0,0.5,0.25. 



dv{t) 
dt 



- -/v{t) + U'ix) =Fo + de{t)v{t) + V2l5CW (21) 



The equation for the energy depot remains instead exactly 
Eg. p5|) . For special case when U'{x) = 0, Fq = and 
e = 0, the corresponding Fokker-Planck equation may be 
solved and the solution is given by Eq. (fTTj) . When includ- 
ing a constant tilt Fq — —a < 0, the stationary Fokker- 
Planck equation may also be solved exactly. The solution 
gives now an asymmetrical distribution with two maxima 
corresponding to the deterministic stable flux velocities 
(see Fig. HI): 



P{i{v) = C exp 



-^--+^log(l + -«^) 



(22) 



The shape of the stationary probability density demon- 



Stationary Probability Density of the Velocity 




Fig. 2. Example of the asymmetric probability distribution of 
the asymptotic velocity under a ratchet in the presence of a tilt 
for different values of the constant force intensity a. The tilt of 
the potential gives rise to the asymmetry of the distribution. 



strated in Fig. 2 shows that under the action of a constant 
force directed to the left (a > 0), flows in both directions 



are possible, however with a different probability: Flows 
directed to the right are less probable than those oriented 
to the left. 



3 Ratchets coupled to energy sources 
3.1 Models of inertia ratchets 

In a series of recent papers [22l[23] various aspects of the 
ABM energetics in an external potential have been an- 
alyzed. In order to investigate further the motion of an 
ensemble of pumped Brownian particles in periodic fields, 
we consider two simple models of smooth 1-dim potentials: 
the symmetric sinusoidal potential (see Fig[3]): 



U{x) = ^ cos(27rx)] 



(23) 



and the asymmetric ratchet potential (See Fig[2]and FigfT^ 
introduced by Mateos and Machura [241125] : 

U{x) = /i{0.499 - 0.453{sin[27r(x + 0.1903)] + 

-|-0.25sin[47r(a;-h 0.1903)]}}. (24) 

We first address the problem of the classical determin- 
istic dynamics of a particle in those potentials. We con- 
struct the bifurcation diagram and identify the origin of 
the current. By analyzing stochastic influences, we detect 
noise-induced directionality of the current. 



3.2 Sinusoidal Ratchets in the trapped regime 

We study first oscillations around x = 0. In this case the 
potential can be approximated by a parabolic one 



U{x) ~ hiT^x^ = — WqX^. 



(25) 



and for a positive value of the bifurcation parameter /5 > 
the system displays self-oscillating solutions. By assuming 
their form as 

v{t) = asin{Lut — a), (26) 
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and substituting into the set of evolution equations, we 
may determine the amplitude a. For small /3 a standard 
derivation of a periodic solution yields a — (2/3)^/^ and 

UJ — UJQ. 

In Fig. 0] the plot of the amplitude of the asymptotic 
velocity a as a function of the parameter q is displayed. 
We see that small amplitudes follow predicted root law 
whereas divergence between analytical and numerical re- 
sults at higher values of q are due to a breakdown of the 
parabolic approximation of the potential. For higher val- 



Amplitude of the Oscillations - Limit Cycle Regime II 



1.4 

12. 
^ 1 
■o 0.8 



E 0.6 



[2(q/Y-c/d)] 
Numerics 




Fig. 4. Amplitude of the limit cycle as a function of driving 
parameter q for different values beyond the critical values (sym- 
bols) in comparison with the theoretical result (full line) .We 
used a ratchet with h — 1, driving parameters 7 — 0.2, c — 0.1, 
d = 0.021. For amplitudes larger than around half of the max- 
imal amplitude, the theory breaks down due to significant de- 
viations of the sine-potential from a parabola. 



ues we observe a bifurcation from a limit cycle to an open 
trajectory. The transition to open trajectories is expected 
in the region where the kinetic energy exceeds the maxi- 
mal potential energy: 



1^2 l(l^c 
2 ° 2\j 



(27) 



Consequently, in what follows, we use the energy trans- 
fer parameter d as an order parameter in our study. The 
equality sign (=) in the above expression defines the value 
of the energy exchange parameter d giving rise to a bifur- 
cation between stable oscillations and the flux regime. We 
denote this value ds- 



dR = 



1 — dch/ c 



(28) 



differs significantly from the value predicted in this ap- 
proximation {dE,a) (See FiglH). 

For particles entrapped in the potential well and per- 
forming sustained oscillations one can derive an analytical 
expression for the value of energy in the energy-reservoir. 
The starting point for the derivation is an assumption that 
at a moderate energy transfer parameterized by d, a con- 
fined closed orbit occurs in a potential well. Let's assume 
that this trajectory is characterized by a sinusoidal veloc- 
ity: 

Xoo{t) = asin{ujt), (29) 

After substituting Xoo{t) in Eq. (|15p and integration, the 
analytical expression for the depot energy e{t) takes on 
the form: 

eoo(t) = (gW^(t)+eo)e-"(*) (30) 



where 



and 



where c?c = jc/l- The above expression, as based on the 
parabolic approximation of the potential, gives only a rough 
estimation of c?^; and is no longer valid in the dynam- 
ical region where the transition from stable oscillations 
(trapped trajectories) to the flux regime takes place. As a 
consequence, the numerically evaluated parameter (d_B.n) 



Wit) 



da^ 



■q{t) = (c + da^/2)t - —sin(2Lut). 



(31) 



(32) 



The integral Ea. ((3T|) can be evaluated in the long time 
limit and the result implemented again to the formula 
Eq. ([50]) leading to an asymptotic expression for the depot 
energy 



q 



da'^/2 



exp 



d^ 



sin(2Lut) 



(33) 



which averaged over the time becomes: 



^00 — ('-oo(^))r 



exp 



sin{2ijjt) ) 

4w 7/, 



c + da^/2 

Here brackets represent the time average 

{f{t))r= lim - r f{t)dt 
T^oc T Jq 

In a first approximatiorQ we can use 



exp I sm(2u}t) 

4cj 



due to the average {sin{2ujt)) r ~ 0. This means that 



da^/2' 



(34) 
(35) 

(36) 
(37) 



^ This approximation holds only for <C 1, which is a 
valid value for our case. A more general value can be found 

as follow: we have that {e^'"'"'^'^^'^^)e[e~^ ,e^], so we can 
put as approximation the mean value between the two extremes 

of the exponential function: (eiEr"'"'^"*^) « cosh{^). In 
this more general case, the relation between the parameters is: 



c+da^/2 



This relation is not simple to invert 



and also requires lo as an input value to obtain a. 
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which differs from Eq. (fT7|) by a factor 1/2. On the other 
hand, by assuming close-to-equilibrium condition (dissi- 
pation term vanishing over the time average) in Eq. (jl4p 
we can evaluate: 



This means that: 



(38) 
(39) 



Comparing Ea. p4)l and Eq. ((39)) . we have a relation be- 
tween the parameters of the equations and the amplitude 
a of the limit velocity of the particle: 



(7 - deoo{t))T = 

- 1 
a 



q 7 ^ /2(? 2c /— 



(40) 



(41) 



This predicted value for a is very close to that found by 
the numerical solution of the equations. From the former 
equation is it possible to evaluate the limit threshold giv- 
ing the stable oscillation of the system. In fact, by putting 
a ~ we obtain the value 

7c 

q ' 

To resume the present result, we have, for dc < d < ds- 

— stable oscillation of v{t) with amplitude a 

— stable oscillation of x{t) in a certain well of the poten- 
tial with amplitude ajlixu 

— oscillating stationary behavior for eoo(i) with frequency 
We = 2a;, 

while for d < dc we have: 

— limit value x(t) in a certain final position Xf 

— decreasing of v{t) {a — 0) 

— damped oscillating behavior with saturation for eoo{t) 
(see Eq.® for = 0) 

Results of numerical evaluation, together with the analyt- 
ical result Eq. p3|) . are plot in two figures below (FigsO 
[6]). As it can be deduced from the simulations, motion in 
the periodic ratchet potential leads to oscillatory behav- 
ior of the energy e{t). That refrains us from using strictly 
the elimination scheme based on the assumption e{t) = 0. 
Instead, the periodic variation of e(t) leads to an aver- 
age (over time) constant value of Coo that differs from the 
value predicted in Eq. (fT7|) . 

3.3 Sinusoidal ratchets in the flux regime 

We will try first a perturbation theory developing the ex- 
pression of the velocity of our system: 



V — Vq + Vi{t) -f 

e = Co + ei(t) 



(42) 

For large driving and small forces (U'{x) « 0) the particles 
move as free and there are two attractors of the velocity 

7 
d 



eo 



(43) 



Xq = 0.0, V(, = 6.0, Co = 2.0 I d = 0.024, y = 0.20, c = 0. 1 

1.5 I 1 1 1 1 1 1 31.2 

Numerical 





9.1 




Numerical 


9.096 






9.092 






9.0SS 






9.084 







2^. Graphs of 



Time Time 

Fig. 5. Asymptotic behavior of the system in the trapped 

case with limit cycle (dc < d < ds)- a 

time-dependent energy represent variation of the depot energy 
e{t). 

Xo = 0.0, Vq = 6.0, eo = 2.0 I d = 0.012, y = 0.20, c = 0. 1 

1 30 




25 
20 
X 15 
10 



10 20 30 40 50 60 70 

Time 




10 20 30 40 50 60 70 80 

Time 



10 20 30 40 50 60 70 80 

Time 



Fig. 6. Asymptotic behavior of the system in the trapped case 
with zero limit (d < dc)and imaginary parameter a. Graphs of 
time-dependent energy represent variation of the depot energy 
e{t). 



We take these solutions as the first term in a perturbation 
series. Inserting this zero- approximation into the Eq. I14[ 
we have: 

. , , dvAx) dx dvi(x) n , . 

«i {t) - — ^ -77 ~ — = -U' (x) (44) 
ax at ax 

^From which, integrating in x, we have the solution: 

Vi{x)vo ~ -U{x) + Const (45) 

and therefore, up to the first order: 

Uix) Const 
^v^ + vi{x)f=iv^) 1 . (46) 



The choice of the constant, is made in order to have a the 
mean value of the potential equal to (See Eq. ([24| ). The 
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value of the constant is then Const — 0.499/i. With this 
choice we obtain, up to the 1st order of approximation: 



= 0.26, y= 0.20, c = 0.1 



U{x) 0.499/1 



Wo 



(47) 



The above expression is valid, in principle, for any shape 
of differentiable potential ratchets. 

Taking into account the Eq. 271 we can make the as- 
sumption for our asymptotic velocity: 



ioo{i) = b + asin{ujt), 
the analytical expression for the e{t) is then: 



where 



and 



with 



,{t) ^ (qZ(i) +eoe-2'iW-')e-«'W 



Z{t) = I e^^'">dt' 
'o 



C{t) = 2dbacos{ujt) 



da" 



-sin[2ujt). 



(48) 

(49) 
(50) 
(51) 
(52) 



Because the term {c+ db"^ + da? / 2) is greater than zero, the 
equation has a non vanishing and non diverging asymp- 
totic behavior. The asymptotic expression we can extract 
from the equation is now: 



eoo(i) 



c + db^ + da?/2 



(53) 



The expression for eoo{t) presents in the exponential two 
oscillating functions with two different frequencies (w and 
2uj). Because the coefficient of the cosine function is greater 
than that of the sine, the observable frequency in the reser- 
voir energy is the same than the velocity one, while in the 
trapped case we have only the frequency 2uj for the asymp- 
totic energy. FiglT] shows the behavior of this expression 
compared with that one obtained by numerical evaluation. 
The agreement is very good. The parameter a and b can be 
evaluated by means of the comparison with the Eqs. (j48| 
and (|47p . while lo has been extracted from the numerical 
evaluation of the limit velocity Voo{t). 




Fig. 7. Case of asymptotical flux. Numerical and analytical 
behaviors are quite superimposed. Differently from the trapped 
case, the energy oscillates mainly with the same frequency of 
the velocity. 



Bifurcation Diagram 
x„= 0.0, v„ = 6.0, e„ = 2.0 I q = 1.0, Y= 0.2, c = 0.1, h = 2.0 
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Fig. 8. Bifurcation diagram for the amplitude (left hand side) 
and the mean value (right) of the asymptotic velocity. The 
comparison between the predicted values and the measured 
ones are plotted. The three region I, II, III, represents the dif- 
ferent dynamical regimes as a function of the order parameter 
d. The indicated threshold values dE,n dE,a represents respec- 
tively the numerical and analytical evaluation (assuming the 
harmonic potential profile) of the escape bifurcation parameter 
ds (see Eq|28ll 



3.4 Influence of external forces and noise 

Assuming that the external force is different from zero, we 
put a certain bias to the right or left direction. A particle 
which is able to go uphill at the cost of the supply from the 
reservoir energy may convert this energy into mechanical 
or electrical energy and perform work. This simple mo- 
tor device is of particular interest in nanobiotechnology, 
where the reservoir energy is in most realistic cases just 
the chemical energy of reaction. 

Among various ratcheting devices, a special class are 
so called staircase ratchets. This type of ratchets may find 



applications to understand and design biological stepmo- 
tors [15,17 . Usually stepmotors have two phases: (i) a 
power stroke where the legs move against a force and (ii) 
a phase of free motion or diffusion. We propose to model 
those by a staircase ratchet having a steep and a flat re- 
gion. 

As a simple example we study the potential obtained by 
a sinusoidal ratchet with a very strong tilt (See the first 
plot in Figs.[9landfT0l): 



U{x) = h 



X sin(27ra:) 

27r ^ ^ 



(54) 
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Here the parameter h denotes the height of one step. As 
can be inferred from Figs. [9] and [TOl the readjustment of 
the exchange energy parameter d can drive the system 
uphill in the staircase ratchet. Fig. [10] presents a longer 
stabilization time because we are close to the critical dcr 
value below which the uphill motion is not observed. In 
fact, for d = 1.3 we observe only a sliding down motion of 
the system. Moreover, in the chosen potential the initial 
velocity plays a key role. Initial velocities lower than the 
limit mean value (equal to 1 in Figs [9] and [T0|) , give rise 
to slopping down motion only. 



h = 1, xj = 0.0, Vo = 2.0, Co = 0.0 I d = 2.2, y= 0.2, c = 0.1 




0.5 1 1.5 2 2.5 3 3.5 4 12 12.5 13 13.5 14 14.5 15 15.5 16 

X X 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 IS 20 



Time Time 

Fig. 9. Active motion on staircase ratchets for a rather large 
value of the strength of driving d = 2.2. 



h= l,x„ = 0.0, Vo = 2.0, eo = 0.0 ld = 1.4, y= 0.2, c = 0.1 




0.5 1 1.5 2 2.5 3 3.5 4 12 12.5 13 13.5 14 14.5 15 15.5 16 

X X 




2 4 6 8 10 12 14 16 18 20 2 4 6 8 10 12 14 16 18 20 

Time Time 

Fig. 10. Active motion on staircase ratchets for a smaller value 
of the driving which is near to the threshold of uphill motion 
d = 1.4. 



the equation motion. As expected, different white noise 
intensities don't lead to any asymmetry in the statistics 
of the flux directionality (right or left) in the absence of 
the constant force {Fq = 0). In contrast, the additional 
constant force, strongly affects the motion. For very low 
values of noise, such as in the deterministic regime, we ob- 
serve that superposition of both forces added gives rise to 
different flux directionality for the set of parameters used. 

When increasing the noise intensity in the region D « 
[10~^, 10~^], a non-monotonic behavior of the fraction of 
right-oriented trajectories Nuight/Nxot as a function of 
the noise intensity is observed. By further increasing the 
noise D « [10~^, 10~^], a region of strong preference to- 
wards left-oriented trajectories is detectable, giving the 
expected behavior of the flux in agreement with the ap- 
plied constant force. Still higher values of noise intensity 
blurred the picture and the rectifying propertie sof the 
system are no longer detectable [N Right /Ntoi ~ 0.5). 

This stochastically driven directionality of the sinu- 
soidal ratchet system in the presence of an additional tilt 
appears to be a generic effect and will be further dis- 
cussed (for asymmetric ratchets) in the forthcoming sec- 
tion, where similar statistics of Njnght/NTot events has 
been evaluated. 



Right Trajectories in Sinusoidal Potential 




lO"'* 10"'^ 10"'° 10* lO"' 10* 



Noise (D) 

Fig. 11. Statistics of the right/left trajectories for the sinu- 
soidal ratchet with and without the presence different constant 
forces as a function of noise-intensity. 



4 Asymmetric Smooth Ratchet Potential 
4.1 Analytical results from perturbation theories 



As in the case of the sinusoidal potential, also for the 
asymmetric ratchet case we observe three different regimes 
As a last result for the sinusoidal- like ratchet, we present depending on the depot-energy influx to the mechanical 
numerically evaluated statistics of the right trajectories as motion: 

a function of the noise intensity for both: the tilted and (i) the state of complete rest (point attractor of the dy- 
tilt-free ratchet. In FiglTT]we plot the case of sinusoidal namics), 

potential with and without additional constant force in (ii) sustained oscillations in one well (bounded attractor) , 
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Ratchet Potential 





1 1 1 


h= i.bo 






h = 0.50 






'^^^^^^^ h = 0.25-"^— — ^ 















Fig. 12. Ratchet potential obtained adding two sinusoidal 
functions for different height h = f.O, 0.5, 0.25 (Eq.[24}. 



(iii) flux regime (open attractors). 

The bifurcation plot of these regimes is isomorphic to that 
one of the sinusoidal potential and it is shown in FiglH We 
study here only the flux regime which is the most interest- 
ing regime for rectifying ratchet systems. The analytical 
expressions for the trapped cases ( corresponding to the 
harmonic approximation of the potential), are identical to 
those reported in the previous chapter. 

Examples of trajectories in the flux regime are given in 
Fig ll3ll4l By inspection of these plots we conclude that 
the change in the right and left asymptotical velocity is 
not monotonically dependent on the energy transfer pa- 
rameter d. 

The analytical evaluation of the velocity and the en- 
ergy in the ratchet case, follows the same scheme used for 
the sinusoidal potential. Using the expression of the veloc- 
ity derived in the perturbation treatment (Eq. (|47p ). with 
the ratchet of the Eq. ([24|) , we can proceed with the eval- 
uation of the analytical asymptotic behavior of the depot 
energy using the expression of the velocity written in the 
following way: 

ioo(i) = ^ + asin(cjt) -I- /sin(2cjt), (55) 

the formal analytical expression for the e(<) is then: 

e^(t) = {qY{t) -f eoe-'^''(/+2a)/^)jg-c(t) 



where 



and 



with 



Y{t) 



at) = (c-f -t-^)t- Ai, 



9{t) — 2dbacos{u!t) — daf sm{LLjt) + 
+ dfhcos{2uot) + ida^ sin(2tjt) + 

-I- \dfasm{?>ijjt) + id/^ sin(4a;i). 
3 8 



(56) 
(57) 
(58) 



Because the term (c-f dfc^-f (ia^/2-|-c?/'^/2) is again greater 
than zero, the equation has a non vanishing and non di- 
verging asymptotic behavior. The asymptotic expression 
we can derive from the equation is then: 



■ + d})^ + da?/2 + dp/2 



(60) 



The expression for eoo{t) contains in the exponential six 
oscillating terms with four different frequencies (w, 2a;, 3a;, 
and 4a;). The relative weight of the higher frequencies are 
low because of the coefficient of these terms and also be- 
cause the value of the coefficient / which is equal, in our 
potential, to one forth the coefScient b. The relevant oscil- 
lating term is then that one with frequency w, the others 
representing less significant harmonics. 
Performing the average of the asymptotic energy (eoo(i))r 
as made with the sinusoidal potential, and equalizing it to 
the value 7/d, we can simplify the expression for the en- 
ergy eoo(i) as: 

eoo(<) = 5e«(*)/-. (61) 



that is the expression used in our predictions. Fig fT51 shows 



0.0, eo = 2.00 I d = 0.260, y = 0.20, c = 0. 1 




(59) 



Fig. 13. Asymptotic behavior of the system in the flux case 
with the usage of the Mateos asymmetric ratchet potential 
plotted in Fig|T2] 



the behavior of this expression compared with that one ob- 
tained by numerical evaluation. Even in this asymmetric 
case, the agreement is extremely good. The parameters a, 
6, / come from the expression of the velocity, and oj has 
been extracted from the numerical evaluation of the lim- 
iting velocity Voo{t)- Otherwise, uj can be also predicted 
by means of the mean flux velocity (t;oo(i))r = b using the 
relation: 

b 2nb 

bT = - = = L = l^a; = 27r6, 62 

where T is the period of the oscillations, i/ the correspond- 
ing frequency and L is the length of the ratchet periodicity. 
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4.2 Results of simulations 

Numerical simulations performed by Tilch et al. (1999) 
have documented that, at least for ratchet models with 
piecewise linear potentials, the onset of a directed net 
current appears in two different directions. Similar behav- 
ior is registered in this study: at sufficiently large values 
of driving, close to the stationary states Vq and Vq , the 
dynamical system possesses open attractors correspond- 
ing to the left or right current states. As we can see in 



Left-Right Trajectories 



Ratio of Right Trajectories 






x„= .0,V(,= .00,e„=2.00ld=.160,-i=.20,c=.10 Xo= .0, v„= .00. e(,=2.00 I d= .150. 1= .20. c= .10 

3 I 1 1 1 1 1 1 1 1 3 

2.5 
2 
1.5 
1 

> 



-0.5 
-1 
-1.5 
-2 

-2 -1.5 -1 -0.5 0.5 1 1.5 2 
X 

Fig. 14. Examples of escape events from the well of the 
ratchet-potential. Relatively small variations of the energy 
transfer parameter d can induce different directionality of flux. 



Fig. mi the ratchet driven by active friction with a strong 
depot-particle coupling (d > d^) possesses 2 momentum- 
dependent attractors. Slight variations of the intensity d 
result in changes of the flux directionality, especially in 
the domain of low noise-strengths (cf. Fig. 15) where the 
dynamics follows the trend observable in a deterministic 
system. The overall pattern, however, does not reveal any 
systematic tendency. In case of stochastic dynamics, two 
asymptotical velocities can be compared with two modal 
values (most probable ones) of the stationary probabil- 
ity density Po{v). Here, instead of examining positions of 
maxima of Po{v), we have analyzed portion of trajectories 





' ' ' ' d=0.bo — 




(1=0.58 




d=0.50 _ 




d=0.48 




d=0.40 







10""' 10" 10° 

Noise (D) 



Fig. 15. Fraction of right-directed trajectories Nmght/NTot 
out of the total number of analyzed cases Ntoi = 5 X 10*. Di- 
rectionality of the asymptotic velocity is displayed as a function 
of the noise intensity. 



whose long time evolution results in a preferred direction 
of motion (Fig. (fTS]) ). Closer inspection of Fig. p5|l allows 
to detect noise-induced changes in the directionality of 
current. For very low noise intensity, a deterministic sce- 
nario prevails and the direction of current is decided by the 
deterministic dynamics (cf. Fig. (|14p ). In contrast, at in- 
creasing values of the noise intensity (D ~ [10~^^, 5-10~^]) 
the system exhibits a tendency to move towards right. At 
still higher values of .D ~ [5 • 10~^, 10~^], the preferred 
asymptotic direction changes to the left. Eventually, for 
values of ~ 10^"^ the motion of the system becomes 
fully delocalized resulting in equal ratios of trajectories 
going to the left and right. 



4.3 Tilted ratchets 

First we discuss the existence of an unidirectional ratchet. 
Studying a Mateos ratchet it is allowed a choice of param- 
eters such that the average of the smaller slope (increasing 
left to right) can still be overcome by the driving mech- 
anism. However the large slope (from right to left) is too 
large to be overcome. In other words, there exists an up- 
hill solution for the smaller slope and no uphill solution for 
the larger slope. This prevents any possibility to go left in 
our case. The ratchet-system is acting then as a unipolar, 
rectifying device. 

Under tilted ratchets we understand ratchets with a 
constant average slope. In other words we have a global 
incline of the ratchet which is due to some constant aver- 
age force. This may model a constant external load against 
which the ratchet has to do work. We mention that several 
authors considered also the case of oscillating tilts. We are 
interested in doing work against a load, therefore consider 
here only constant external slopes. 

In Fig. [in] we give an example of an uphill motion 
against a force (here F = 0.03). The load force is in the 
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Trajectories F = +0.03 
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Fig. 16. Examples of uphill flux on tilted ratchets. The load 
force is directed left to right. 



example directed left to right. We see that in dependence 
on the initial conditions both directions may be realized. 



5 Applications 

The investigated systems may be of interest for the mod- 
elling of molecular motors which are able to convert chemi- 
cal energy into mechanical energy (work). We see possible 
applications to proton pmiips, the rotating motors con- 
necting with the work of ATP-ase [T51[n] and also to step 
motors as proposed by Bier [T5jll7j. 

Molecular motors are nanotechnological objects be- 
ing the result of biological evolution. All molecular mo- 
tors use the energy quantum connected with the syn- 
thesis or hydrolysis of the nucleotides ATP/ADP or the 
difference of the electrochemical potentials on the cyto- 
plasm membrane. Standard models of molecular motors 
are based on the Smoluchowski equations for discrete sys- 
tems having several states which correspond to attach- 
ment or detachment p!^ll4j . Many models have been de- 
veloped which follow similar lines. We follow in this work 
another route which is based on Hamiltonian ratchets. We 
studied Hamiltonian ratchets which are connected to an 
energy reservoir and gave special attention to possible ap- 
plications to molecular motors. We investigated in detail 
the motion of a particle against a gradient of the poten- 
tial i.e. uphill motion under conditions where the external 
force is pointing downhill. The general schema is the fol- 
lowing: chemical energy is absorbed in the form of ATP 
and introduced into our "motor" increasing the reservoir 
of e{t) by a certain amount. This is modelled here by a 
continuous inflow q. In some other work we develop a more 
refined model based on the assumption of discrete energy 
quanta, representing the absorption of one molecule ATP 
[18]. The absorbed energy flows to the "motor" and is 
transformed into mechanical or electrical potential energy. 
This could e.g. model the increase of the energy of protons 
by transport through the membrane. 



In the case of the ATP-ase motor, in one of the di- 
rection of its rotational kinematics, the enzyme ATP-ase 
hydrolyzes ATP into ADP and anorganic phosphate-Pi, 
releasing energy and moving protons. The work of FqFi- 
ATP-ase is connected with the rotation of a "rotor". The 
idea about a " rotor" corresponds to our knowledge about 
the structure (see Fig. [T7|l . In our model of active Brown- 




Fig. 17. Scheme of the molecular motor FiFo-ATPase- 
synthase. The complex Fo which is imbedded into the 
membrane is responsible for the transport of protons through 
the membrane. The complex Fi sitting on the membrane is 
responsible for the catalytic functions during the synthesis 
/ hydrolysis of ATP. The figure has been adapted from an 
online publication of the Royal Swedish Academy of Sciences 
(http://nobelprize.org/nobeLprizes/chemistry/laureates/1997/ 
press.html). 



ian particles the variable x has to be interpreted then as 
the rotation angle of the rotor. 

Another possible application may be to model step mo- 
tors, as studied e.g. by Bier [l5lfT6l[T7] . Here the motion 
of molecular legs is composed by two phases. The first 
phase is the power stroke where the system works against 
a force, and the second phase is the free or diffusive mo- 
tion. One possibility to model such two-phase motions is 
to use staircase-like ratchets defined by Eq.(52). Another 
possibility is to develop ratchet models based on two ac- 
tive particles coupled by a spring. Step motor models of 
this type, based on overdamped ratchet motion, were first 
studied by Derenyi and Vicsek pS] . 



12 



Alessandro Fiasconaro et al.: Active Brownian Motion Models and Applications to Ratchets 



13. F. Jiilicher, J. Prost, Phys. Rev. Lett. 75 (1995) 2618, 78 
(1997) 4510. 

14. F. Jiilicher, A. Ajdari, J. Prost, Rev. Mod. Phys. 69 (1997) 
1269. 

15. M. Bier, Phys. Rev. Lett. 91 (2003) 148104-1 

16. M. Bier, Contcrnp. Phys. 46 (2005) 41-51 

17. M. Bier, BioSystem 88 (2007) 301-307 

18. W.Ebeling, E. Gudowska-Nowak, A. Fiasconaro, Acta 
Phys. Pol. B 39 n.5 (2008) 1225. 

19. F. Schweitzer, W. Ebeling, B. Tilch. Phys. Rev. E 64 
(2001) 02110-1-12. 

20. U. Erdmann, W.Ebeling, and V. S. Anishchenko. Phys. 
Rev. E 65 (2002) 061106. 

21. U. Erdmann, W.Ebeling, A. Mikhailov, Phys. Rev. E 71 
(2005) 051904. 

22. B. Tilch, F. Schweitzer, W. Ebeling, Phys. A 273 (1999) 

293-314 

23. F. Schweitzer, B. Tilch, W. Ebeling, Eur. Phys. J. B 14 
(2000) 157-168 

24. J.L. Mateos, Phys. Rev. Lett. 84, (2000) 258; Acta Phys. 
Pol. B 32 (2001) 307-320 

25. L. Machura, J. Luczka, P. Talknor, P. Hariggi, Acta Phys. 
Pol. B 38 (2007) 1855. 

26. A. Derenyi, T. Vicsek, Proc. Natl. Acad. Sci USA 93 
(1996) 6775-6779. 

molecular motors connected to the ATP synthesis/hydrolysis 
have been briefly discussed. 



This work has been supported by the Marie Curie 
TOK grants under the COCOS project (6th EU Frame- 
work Programme, contract No: MTKD-CT-2004-517186) 
as well as the ESF programme STOCHDYN. The authors 
thank Martin Bier for discussions and advice. 



References 

1. M. Smoluchowski, Phys. Z. 17 (1916) 557-571. 

2. Yuri L. Klimontovich, Statistical theory of open systems. 
Kluwer Academic Publ., Amsterdam, 1997. 

3. W. Ebeling, 1. Sokolov, Statistical thermodynamics and 
stochastic systems World Scientific Publ., Singapore 2005. 

4. F. Schweitzer, W. Ebeling, B. Tilch, Phys. Rev. Lett. 80 
(1998) 5044-5047. 

5. W. EbeUng, F. Schweitzer, B. Tilch, BioSystems 49 (1999) 
17-29. 

6. U. Erdmann, W. Ebeling, L. Schimansky-Geier, F. 
Schweitzer, Eur. Phys. J. 15 (2000) 105-113. 

7. T. Vicsek, Fluctuations and scaling in biology. University 
Press, Oxford, 2001. 

8. Alexander S. Mikhailov and Vera Calenbuhr. From Cells to 
Societies. Springer Series in Synergetics. Springer, Berlin, 
2002. 

9. Frank Schweitzer. Brownian Agents and Active Particles. 
Springer Series in Synergetics. Springer, Berlin, 2003. 

10. A. Ordemann, F. Moss, G. Balaszi. Motions of daphnia in 
a light field: random walks with a zooplankton in F. Beck, 
M.T. Hiitt, U. Liittgc (cds.), Springer, 2003. 

11. W. Ebeling, Acta Phys. Pol. B 38 (2007) 1657. 

12. W. Ebeling, L. Schimansky-Geier, Eur. J. Phys. J, in press. 



6 Summary 

In this work we analyzed a mechanical system with in- 
ertia subjected to a dissipative forcing having two terms: 
one passive damping and an 'active' contribution provid- 
ing energy to the system by means of an energy supply 
described in an additional coupled equation. The system 
has been studied under the influence of smooth ratchet 
potentials: symmetrical (sinusoidal potential) and asym- 
metrical (Mateos- type potential and tilted ones). The ana- 
lytical and numerical evaluation of the equations has been 
found for both the asymptotical velocity of the mechan- 
ical system and the depot energy. The system presents 
bifurcations of the asymptotic velocity as a function of 
the energy transfer parameter d. The three regimes found 
correspond to: 1) relaxation in a potential minimum (van- 
ishing motion); 2) oscillating motion in a well (limit cycle), 
and 3) flux motion with two values of the asymptotic ve- 
locity. The motion in the flux regime is then possible in 
two directions, even in the presence of a tilt in the poten- 
tial. The numerical simulations of the system under the 
action of white Gaussian fluctuations show the effect of 
noise-controlled directionality of the motion. 

Possible applications of the ABM system in modelling 



